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Abstract 

We experimentally and theoretically investigate electromagnetic cavity modes in ultra-narrow 
Al-AlO^-Al and Nb-AlO^-Nb long Josephson junctions. Experiments show that the voltage spacing 
between the Fiske steps on the cur rent- volt age characteristics of sub-/im wide and several hundred 
/im long Al-AlO^-Al and Nb-AlO^-Nb Josephson junctions increases when decreasing the width of 
a junction. This effect is explained by stray magnetic fields, which become important for narrow 
junctions. Theoretical estimates of the Fiske step voltage based on a nonlocal wave propagation 
equation are in good agreement with our experimental data. Using the nonlocal model, we de- 
termine the size and mass of a Josephson vortex by means of a variational approach, and relate 
vortex size to the experimentally measured critical magnetic field of the junction. 



PACS numbers: 74.50. +r; 85.25.Cp; 75.70.-i 




FIG. 1: Schematic view of an edge-type long Josephson junction of width w and length L. The 
lines show the magnetic field B, due to vortex, which penetrates the electrodes and extends into 
free space. The extended stray magnetic field creates the additional surface current J5. 

I. INTRODUCTION 



The electromagnetic properties of long Josephson tunnel junctions have been the subject 
of intensive studies over the past four decade a^'^'?'^ . Most investigations of junction dynam- 
ics have been based on local theory. This theory assumes that electromagnetic fields are 
concentrated in the junction barrier, and neglects the distribution of the magnetic and elec- 
tric fields outside the junction. However, in reality the fields penetrate the superconducting 
electrodes and extend into free space (Fig. Local theory breaks down when the fields 
outside the junction contribute significantly to the junction energy. 

The magnetic field inside the superconducting electrodes of a Josephson junction be- 
comes important when its characteristic decay length is of the order of (or larger then) the 
characteristic spatial scale of variation of the Josephson phase if. This type of nonlocality 
we call internal, because it depends only on the field distribution inside the junction and its 
electrodes, and neglects extended, stray, fields outside the junction. 

When the energy of stray fields becomes comparable to the energy of the field inside the 
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junction, one has to take stray fields into account. These stray fields generate an additional 
surface current J5, which can be represented as a functional of the magnetic field distribution 
inside the junction, leading to a nonlocal equation. Because these effects depend on the 
geometry of the junction, we call them geometrical nonlocal effects. 

Detailed knowledge of Josephson vortex dynamics is important for understanding the 
magnetic flux motion and related phenomena in superconductors. In a long junction, a 
vortex behaves like a classical particle, which is characterized by the spatial coordinate 
of its center and the effective mass. The interest to experimentally verify theoretically 
predicted nonlocal corrections to the conventional model has been biased by recent studies 
of long ultra-narrow junctions where Josephson vortices behave as macroscopic quantum 
objects^. In very narrow junctions, the effective dynamical mass of the vortex is expected 
to be influenced by nonlocally described stray magnetic flelds. Most of the theoretical papers 
assume, however, that the effective dynamical mass of a Josephson vortex is proportional to 
the width of the junction as expected in the local theor y '''1^1^ . 

In this paper, we experimentally verify that, in very narrow long Josephson junctions, one 
has to indeed take into account the vortex mass correction due to the stray flelds outside the 
junction. Signiflcant deviations from the local theory have previously been observed in the 
current-voltage characteristics of long junctions of small width^-. We systematically study 
these effects for different junction materials, and compare the experimental flndings with 
the nonlocal sine-Gordon model. 

A Josephson junction can be considered quasi-one-dimensional, as long as its transverse 
dimension, the width w in z-direction (see Fig.Q), is smaller than the Josephson penetration 
depth 



where $0 is the magnetic flux quantum, is the critical current density across the tunnel 
barrier, A = 2Xi + d is the magnetic thickness of the junction, is the London penetration 
depth and d is the barrier thickness (typically, d <^ X^). In this limit {w <^ Aj), the 
long junction is described by the perturbed one-dimensional sine-Gordon equation for the 
Josephson phase difference ip 
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where subscripts mean partial derivatives with respect to the spatial coordinate x along the 
junction and time t. The parameter 



is the Josephson plasma frequency, C is the specific capacitance of the junction per unit area, 
and a is the dissipation coefficient due to the quasiparticle tunnelling across the barrier. The 
bias current density 7 = j/jc is normalized to the critical current density jc and is, in general, 
dependent on x. 

Equation is derived using the following assumptions: i) ^ w ^ Aj and ii) hu ^ A, 
where 00 is the oscillation frequency of the phase if, and A is the superconducting energy 
gap of the electrodes. The breakdown of any of these conditions requires modification of 
Eq. (j2I). As mentioned above, for large London penetration depths A^, > Aj (or small 
Josephson penetration depths Aj) the magnetic field inside the superconductor starts to 
play an important role. On the other hand, for sufficiently small junction width w ~ Al, the 
stray fields outside the junction have to be taken into account. In either case, the dynamics 
of the phase ip has to be described by an integro-differential equation, i.e. the problem 
becomes nonlocal. When the oscillation frequency u is comparable to the superconducting 
gap A, the frequency dependence of the London penetration depth Xl{^^) according to the 
microscopic theory {material dispersion) should also be taken into accounlpJ-i^iii. 

The nonlocal equation describing Josephson phase dynamics can be written in the general 
form 



where the function Q{x, x') is the nonlocality kernel, which needs to be specifically deter- 
mined for each nonlocal problem. In the local case Q{x, x') = 6{x — x'), where 6{x) is a delta 
function. Several theoretical approaches have been proposed for different types of nonlocal- 
ity, in the case of an infinitely long Josephson junction. In the limit of an infinitely long 
junction, the kernel of the integral in Eq. is reduced to the function of a single argument. 



Nonlocal models can be divided into two groups, those treating internal nonlocality inside 
bulk junctionsi^ii^, and those dealing with nonlocal effects due to outer stray fields resulting 



from the geometry of the junction and its electrode o^^i'^'''i'^^i'^^i^°i^'^i^^i^^i^'^i^^i^^i^'''i^^i^^ . The latter 



group of theories takes into account the field configuration not only inside the junction. 




(3) 





Q{x,x') = Q{\x-x'\). 
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but also around it, and incorporates the finite size of tlie sample and its shape-dependent 
magnetic properties. The shorter the junction is in x-direction (see Fig.Q]), or the narrower it 
is in z- direction, the larger the effect of the geometrical nonlocality. Four typical geometries 
of long Josephson junctions are presented in Fig. El These are (a) edge-type, (b) ramp-type, 
(c) window and (d) overlap long junctions. In the ideal case of an infinitely long junction, the 
electro dynamic problem posed by (a)-(c) type junctions have been solved. In particular, 
theoretical models for the geometrical nonlocality have analyzed the following cases: 1) 
edge-type junctions between thick films w > Xl, neglecting internal nonlocal effects^iiSS; 2) 
edge-type junctions between thin films w < Xl (Pearl's limit 3) edge-type and 4) 
ramp type junctions between films of arbitrary width, taking into account both the internal 
and external nonlocal problems^i*^; 5) a variable thickness bridge above a ground plan o^^i^^ 
and 6) window junctions^SiSLSI. 




FIG. 2: Different geometries of long Josephson junctions formed between two superconducting 
films: (a) edge-type junction; (b) ramp-type junction; (c) window junction; (d) overlap junction. 

In this paper, we study the effect of geometrical nonlocality by analyzing experiments 
on extremely narrow long Josephson junctions, down to a width w of 0.1 fim. The paper 
is organized as follows. In Sec. |n] we describe the geometries and preparation process of 
ultra-narrow low-T^ long Josephson junctions. In Sec. II I II we present experimental results 
obtained for Nb-AlO^^-Nb and Al-A10a;-Al junctions. Theoretical models based on nonlocal 
electrodynamics, and their comparison with experiment, are presented in Sec. IIVI Section fVl 
contains concluding remarks. 
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II. FABRICATION AND GEOMETRY OF ULTRA-NARROW LONG JOSEPH- 
SON JUNCTIONS 



The fabrication process of high quahty long Josephson junctions in order to study the 
effects of nonlocal electrodynamics should satisfy the following three requirements: (i) ability 
to vary the junction width, from several micrometers down to several hundred nanometers, 
(ii) constant width along the junction length L ^ Aj, which is typically several hundred 
micrometers, and (iii) an idle area of overlapping electrodes (window) that is as small as 
possible, so as to prevent it influencing wave propagation in the junction. We prepared 
samples using two different technologies based on aluminium and niobium. All samples 
were fabricated on a thermally oxidized Si substrates. 



A. Al-AlO^-Al junctions 

The Al-AlO^-Al junctions were prepared with the shadow evaporation technique. This 
method for the preparation of sub-/im A1-A10^.-A1 tunnel junctions is well known (see, e.g. 



Ref. l3ll). However, we cannot apply this method directly. A hanging bridge of electron resist 
(e.g. PMMA) or another material is usually used for shadowing during the evaporation of 
Al. As the length of the bridge becomes longer the bridge starts to sag. Thus, for long 
junctions (several hundred micrometers) the bridge cannot be well fixed. Therefore, in our 
preparation method, we use a "shadowing window". Schematically the fabrication steps are 
shown in Fig. |21 

The first step was the formation of Au electrodes, used for spatially-uniform bias current 
injection and voltage measurement leads (Fig. IHla). These electrodes were formed with the 
help of electron-beam lithography, thermal evaporation and the lift-off technique. The Cr 
layer under the Au layer improves Au adhesion to the substrate. The thicknesses of Cr and 
Au layers were 10 nm and 60 nm, respectively. The electrodes for current injection and 
voltage measurement across the junction have different geometries. One sees in Fig. |3] that 
the bias current leads have a length almost equal to that of the Josephson junction. The 
potential leads are much smaller, and are connected to the junction only at the edges of the 
Al electrodes. This geometry provides homogeneous current injection into long Josephson 
junctions and allows four point measurements of junction characteristics. 
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FIG. 3: Schematic diagram of the fabrication procedure for Al junctions, (a) Deposition of Au 
electrodes; (b) formation of a window in the PMMA resist for shadow evaporation; (c,d) deposition 
of top and bottom Al electrodes. 



To form long Josephson junction by shadow evaporation, we used double layer resist 
PMMA-MA/PMMA 950 K, which was spanned onto the substrate with prepared electrodes. 
Then, using electron-beam lithography, we opened a window between the Au electrodes, as 
shown in Fig. Elb. In the next step, we evaporated Al onto a tilted substrate (see Fig. |Hlc). 
The thickness of this bottom Al layer was about 50 nm. Immediately after Al evaporation 
O2 was injected into the vacuum chamber. The pressure of O2 was kept at 10~^ mBar for 
6 minutes. Then, a 200 nm thick top layer of Al was evaporated onto the tilted substrate 
(Fig. Eld). The overlap of the bottom and top Al electrodes defines the junction area. A 
scanning electron micrograph of the resulting Josephson junction and electrodes is shown in 
Fig.ia 

Schematically, the cross section of the fabricated Josephson junction is shown in Fig. Ela. 
Note that the total junction width is composed of two sections, one, of width w', parallel 
to the film surface plane and the other, of width w", which is tilted. The total width is 
therefore w = w' + w" . Typically the relation w" <^ w' is well fulfilled. However, in the 
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FIG. 4: SEM picture of Al-AlO^^^-Al junction. The inset shows an enlarged view of the connection 
of current and potential electrodes to the junction. 
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FIG. 5: Schematic cross sections of experimentally studied Al-AlO^^-Al (a) and Nb-A10x--Nb (b) 
junctions. 

case of the most narrow junctions, the contribution of w" to the junction width becomes 
significant. We estimate w" ~ 60 nm. 



B. Nb-A10a;-Nb junctions 

There are several different approaches for the preparation of Josephson junctions using 
Nb-A10a;-Nb trilayer technology. They differ mainly in patterning and insulation of the 
junction. Junctions with dimensions less than 1 /im usually have a window geometry (see 
Fig.|21c), and are produced by reactive ion etching and vapor deposition of a dielectric layer. 
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^^^^ substrate [^^^^^^^^ cross-linked PMMA 
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\//////A topNb 

FIG. 6: Schematic diagram of the fabrication procedure for Nb junctions, (a) Formation of the 
trilayer; (b) insulation of the trilayer edge; (c) deposition of Nb wiring; (d) removal of Nb wiring 
and top Nb from all areas except that of the junction, (e,f ) narrowing of the junction using electron 
beam lithography and reactive ion etching (repeated several times). 

Using this approach, one could prepare sub /im junctions, however they would not satisfy 
the requirement of small idle region between the superconducting electrodes around the 
junction. In our group, we developed an alternative method, which allows junctions on the 
very edge of the Nb-AlO^-Nb trilayei^. Schematically, this process is shown in Fig. IHl The 
junction cross section prepared in this way is shown in Fig. Elb. 

In the first step, the Nb-AlO^^-Nb trilayer region is defined by lithography and reactive ion 
etching (Fig.lHla). The next step is the fine definition of one of the edges of the structure using 
electron beam lithography and etching. This particular edge was insulated by cross-linked 
PMMA (see Fig. Elb). The details of the cross- linking procedure and the self-alignment of 
cross-linked PMMA along the trilayer edge can be found in Ref. Q. In the next step, the 
Nb wiring layer was deposited as shown in Fig. IHlc. Electron-beam lithography and reactive 
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TABLE I: Parameters of the measured junctions. 



junction material 


Al, nm 


Jc, A/cm^ 


A J, fim 


c, m/s, 


Fgap, mV 


Al (T = 0.3 K) 


15 


270 


56 


8-10*^ 


0.35 


Al (T = 0.9 K) 


17 


160 


68 




0.29 


Nb (T = 4.2 K) 


90 


210 


25 


7-106 


2.7 



ion etching were used to define the junction area (see Fig. IHld). After performing a full set 
of measurements on the junction, the width of the junction was iteratively decreased using 
electron beam lithography and reactive ion etching, as shown in Fig. IHl e and f . 

III. EXPERIMENT 

In this section, we present the data acquired for Al-A10a.-Al and Nb-A10a;-Nb ultra- 
narrow long Josephson junctions. A magnetic field H of up to 40 Oe was applied in the 
plane of the junction, perpendicular to its long dimension. The magnetic field of the Earth 
was screened out with a cryoperm shield. The experimentally determined parameters of the 
junctions are summarized in Table H] The London penetration depths for Al and Nb were 
found from period AH of Ic{H) at high fields for widest junctions, as Al = (^q/ {2LAH). 
The Josephson penetration depth was found from the critical current density using Eq. (P). 
The Swihart velocity was found by fitting Fiske step voltages to the theoretical dependence 
(see Sec. IIVC|1 . The junction width was measured using a scanning electron microscope. 

In the following, we describe two main effects observed in our ultra-narrow Josephson 
junctions, when the junction width w is decreased. These are 1) a reduction in the first 
critical field Hd, and 2) an increase of the Fiske step voltage spacing AVps- According 
to standard (local) theory, the first critical field ifci and the Fiske step voltage should not 
depend on the junction width^. 

The first critical field corresponds to the field when a single vortex enters the junction. 
In the local theory, the first critical field Hd is inversely proportional to the Josephson 
penetration depth Aj. In the following sections we present the experimental data for Hd in 
terms of an effective length scale 

= vri/.i(2AL + t)- 
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Note that this equation is similar to the relation between Aj and Hd in the local theory. 

Fiske steps arise when the Josephson oscillation frequency resonates with cavity modes 
of the junction. The voltages of these steps according to local theory are given by 



where c is the Swihart velocity and n is the step number. Note that the Fiske step voltages 
Vn are also determined by the length scale Aj, provided Josephson plasma frequency does 
not change. 

A. Al-AlO^-Al junctions 

We measured a series of six 230 /im long Al-AlOx-Al junctions. The width w of the 
junctions varied between 0.1 /im and 1 /im. The homogeneity of the junctions was verified 
by measuring the modulation pattern of the critical current Ic vs. field H aX T = 0.9 K. In 
Fig. [71 Ic{H) patterns of three samples of width w = 0.11, 0.25 and 0.82 /im are presented. 
We find the first critical field by linearly extrapolating the central peak in the critical current 
modulation pattern (see inset of Fig.Ej). We observe that the first critical field Hd decreases 
with decreasing the junction width w. 

The experimentally determined effective length scale A j is presented in Fig. |H1 as black 
squares. This length scale depends strongly on the junction width. Theoretical curves are 
discussed below, in Sec. IIVDI 

Fiske steps of all junctions were measured by tracing the current-voltage characteristics 
during a continuous sweep of the external field H in the range from —40 Oe to +40 Oe. 
These measurements were performed at temperature T = 0.3 K, because at the higher 
temperature T = 0.9 K high damping suppresses the Fiske resonances. At T = 0.3 K the 
first critical field Hd of the wide junctions {w ~ 0.6 — 0.8 /xm) was larger than the maximum 
field we could apply. In Fig. IHl the superimposed characteristics are plotted for the junction 
of width w = 0.82 /im. 

We have observed that the positions of the Fiske steps strongly depend on the width of 
the junction. This dependence for the first step is shown in Fig. ^Jand compared with the 
theory presented below in Sec. IIV CI 
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FIG. 7: Critical current modulation Ic{H) in magnetic field for long Al-AlO^-Al Josephson junc- 
tions at T = 0.9 K. The inset zooms in on the modulation patterns of the junctions of width 
w = 0.11 ;um and w = 0.25 /zm. 

B. Nb-A10a;-Nb junctions 

Results of systematic measurements of ten long narrow Nb-AlO^-Nb junctions were pre- 
sented earlier in Ref. 3 We have now additionally measured three more junctions. Mea- 
surements were performed at temperature T = 4.2 K. In the present paper we analyze earlier 
data and our new results more extensively. The junction length is 200 /xm and the width 
w ranges from 0.3 /im to 4.5 /im. The Josephson penetration depth was estimated from 
the critical current density using Eq. (P) to be about Aj ~ 25 yum. The critical current vs. 
magnetic field patterns, reported in Ref. ^ and those of the new junctions show that all 
junctions are rather homogeneous. As in the case of the Al-AlOx-Al junctions, a decrease 
in the first critical field Hd (Fig. ITT|l . and an increase of the Fiske step voltage spacing 
AVps (Fig- EH), with decreasing width w of the junction were observed in the Nb-A10a;-Nb 
junctions. 
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FIG. 8: Characteristic spatial scale of the Josephson phase in Al-AlO^^-Al junctions. Squares 
indicate the characteristic length scale Aj, as extracted from the measured critical field Hd at 
temperature T = 0.9 K, using to Eq. Q. Lines show theoretical estimations of the size Aj of a 
vortex, for details see Sec. IIVDI 

IV. COMPARISON WITH THEORY 



In this section, we present a short review of the existing theoretical models for the nonlocal 
effects in edge-type Josephson junctions (Fig. [T]). In particular, we are interested in the 
nonlocal kernels Q{x,x') to be used in Eq. and in the dispersion relation uj{k) of the 
junction, which determines the voltage position Vn of Fiske resonances. To compare with 
the experimental data, we include the frequency dependence of Xl^^^) (material dispersion), 
and relate the critical magnetic field Hd to the effective size of a vortex in the quasi-local 
limit. 



A. Geometrical nonlocality 

In a junction of finite length L the translational invariance along the junction is broken 
and, in general, the kernel Q{x, x') depends not only on a; — x', but also on the sum x + x' 
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FIG. 9: Superimposed traces of Fiske resonarices of the current-voltcige chctracteristics of Al-AlO^;- 
Al long Josephson junctions at several different magnetic fields. The junction width is = 0.82 ^m 
and the temperature is T = 0.3 K. 

and the length L. However, in a finite junction, deviations of the nonlocahty kernel from 
the one of an infinitely long junction are expected, mainly in the edge region on the order 
of the magnetic screening length Al (or \p = 2\\/w in thin films). These deviations can 
be neglected for the long junctions L ^ \l {L ^ \p) considered in the following. This 
assumption is justified by the fairly good agreement obtained with the experimental data 
presented in figures H CUl HH [H [121 and d 

We distinguish between nonlocal effects due to stray fields (and screening currents) inside 
and outside the junction electrodes. These two contributions are referred as internal and 
geometrical nonlocality, respectively. 

(a) Let us first consider the theory for internal nonlocal effects in bulk iunctions of width 

nn 

w ^ \l, see. Refs. I14J)15L In this case the kernel in Eq. (jH) is 



Q{X) = {l/Ti\L)M\x\/\L), 



(7) 



where Kq is a modified Bessel function, having a logarithmic pole at x — and exponentially 
decaying at large distances^Si. Therefore, nonlocal effects become important here, when 
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Junction width (^m) 

FIG. 10: The width dependence of the first Fiske step voltage Vi for Al-AlO^-AI junctions. Symbols 
correspond to experimental data obtained at T = 0.3K. The solid line corresponds to the theory 
of case2i (c) and the dashed line to case^^ (d). 

is large compared to the typical length scale l^, on which the phase varies , i.e. in the case 
of Xl > Aj. In the limit \l ^ Aj the local sine-Gordon equation is recovered. 

Usually A J is much larger than A^, because the critical current density jc in the junction 
is much smaller than the critical (depairing) current density jd = 2$o/ (3-\/37r^/ioA|^) in the 
bulk superconductor is the coherence length). Exceptions are, for example, junctions in 
high-temperature superconductors created by planar defects such as twins, stacking faults, 
low-angle grain boundaries, eto^. These structural defects often form Josephson junctions 
with high jc and therefore small Aj. In terms of the critical current density across the 
junction, the condition Aj < A^ is fulfilled if jd/ k < jc < jd, where k = Xl/C, is the 
Ginzburg-Landau parameter-^. Note that for extreme type-II superconductors, k ^ 1, the 
relation Aj < A^ holds over a wide range of jc- 

Linearizing the wave equation {simp ~ one obtains the dispersion relation uj{k) for 
small-amplitude linear electromagnetic waves (p{x,t) = ipoexp{—iLjt + ikx) (here \ipo\ <^ 1) 
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FIG. 11: Characteristic spatial scale of the Josephson phase in Nb-AlO^-Nb junction at T = 4.2 
K. The open circles denote the length scale Aj extracted from the experimental critical field Hd 
according to Eq. © for new junctions. The solid symbols correspond to experimental data for Hd 
of R.ef.llCI. Lines correspond to the theoretical estimations of the size Aj of a vortex (see Sec lIVDj) . 



propagating along a Josephson junction. From the above kernel it follows that 



1 + 



(8) 



(b) The geometrical nonlocality of an edge-type junction formed between two thin su- 
perconducting films w < Xl (Fig. |2l(a)) has been extensively studied^'^^iSS. In this case, the 
stray fields outside the junction area account for the entire electromagnetic energy of the 
junction. The typical length scale upon which magnetic fields vary, is not the London pen- 
etration depth Xl, but the Pearl penetration depth^P Xp = 2X\/w ^ Xl. The nonlocality 
kernel is then given by2& 

dkdq exp(iA;x) 



Q{x) 



4Az 
w 
1 

2Xt 



(27r)2 ^WT^[l + A|,g21 



Hn I ^ 

Xp 



(9) 



where Hq and Yq are the Struve and Bessel function of the second kindfl^". As in the previous 
case, this kernel has a logarithmic pole at x —>■ 0, however at large distances it decays as 
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FIG. 12: The width dependence of the maximum of the first Fiske step voltage Vi for Nb-AlO^^-Nb 
junctions. The open symbol correspond to new experimental data and the data for solid symbols 
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are taken from Ref. llOl . The solid line is calculated using the theory of Ref. 



corresponds to the theory of Ref. 
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and the dashed line 



~ 1/x. When the phase variation length is large compared to the magnetic screening 
length, Ap ^ /<^, the characteristic vortex size is given by \j = = Xj^Xl/wY^"^ ^ Xj. 
The limit of Ap ^ l^^ is then equivalent to Xl ^ Xj{w/ XlY^"^ . In the opposite limit of 
lyy -C Ap or Xl ^ Xj{w/ XlY^'^ the length scales for the tails (Aj) and the core (Aj/A^) of a 
Josephson vortex are differenli^. 

(c) Ivanchenko's theorj*2i applies to edge-type junctions of width larger than the London 
penetration depth, w > Xl. This model assumes that stray fields outside the film affect only 
the surface, and that the interior of the junction is described by the local theory. In this 
case the integral kernel has the form^i: 

AXl f dkdq exp(ifcx) 



Q{x) = 6{x) + 



w 



(27r)2 ^/pT^[1 + Aig2] 



*(^) - ^ (h» (it) - (it 



(10) 



Note that the second term is identical to the kernel of Eq. 0, corresponding to the case of a 



17 



thin film junction, but the Pearl penetration depth Ap is replaced by the London penetration 
depth A^. As in the previous case, when l^p is large compared to the typical length scale of 
the kernel, A^ ^ l^, the vortex size is given by Aj = = Aj(l + Xl/wY^"^- 

In the limit of ^ Al for cases (b) and (c), Eq. does not transform exactly into the 
local sine-Gordon equation with effective length scale Aj, because Q{x) does not correspond 
to 5{x) for Xp/L — > 0. Nevertheless, the locality of the kernel suggests that in this limit 
an effective local sine-Gordon equation with renormalized Josephson penetration depth Xj 
should be qualitatively correct. Then one could use the expressions for Hd and c = cUpAj 
from local theory, if Aj is replaced by a new length Aj, see below. 

In the absence of dissipation and bias current, using the kernel we obtain the dis- 
persion relation for this case as 




cj(fc) = cjp./l + PA^ + ^-^^ arccos ( ) (11) 



for kK > 2 and 



for kA < 2. 

(d) The theory developed by Lomtev and Kuzovlev^^ is the most general as it contains 
the special cases (a)-(c) discussed above. It assumes that the phase ip{x) does not vary over 
the width w. The kernel is 

irXLyxJ^'KwXlJo K^{K + qcoth{Kw/2)y ^ ^ 

where k = + A|)^/^, Jq is the Bessel function of the first kind. In this case the dispersion 
relation in the absence of perturbation terms in Eq. ^ is24 



/ U2\2 1.2 \2 

u{k) = u 1 + ^fJ^ + ^Fik) , (14) 

y a/1 + FAi ttAl 

where 

oo 

2 f 1 1 1 

^^^^ K + k coth{Kw/2) (g2 - k^y/2 ■ 

k 

This theory has also been generalized to describe ramp-type junctions^^ (Fig. |21^b)). 



18 



Let us quote typical parameters for our low-Tc long Josephson junctions. For the Nb- 
AlO^^-Nb junctions considered in this paper Aj ^ 25 fim, and ^ 90 nm at T = 4.2 K, 
and for Al-AlO^^-Al junctions Aj ~ 56 /xm and A/, ~ 15 nm at T=0.3 K, and Aj ~ 68 fim 
and Al ~ 17 nm at T=0.9 K. For these parameters, the condition Al <^ Aj is fulfilled in 
all cases, i.e. the internal nonlocal effects as discussed in case (a) are negligibly small. For 
Al-A10a;-Al junctions w ^ 0.1 — 0.8 fim, and for Nb-AlO^^-Nb junctions w ~ 0.3 — 4.2 fim, 
which is significantly larger than Xl. It is reasonable to expect that our junctions can be 



approximately described by the theory presented in Ref. |21|. In the following analysis, we 
compare our experimental data with Ivanchenko's theory^, case (c) and with the more 
general theory of Lomtev and Kuzovle'^^, case (d). 



B. Material dispersion 

We will see in Sec. II V C\ that for small frequencies the nonlocal dispersion relation 
predicts a strong decrease of the spacing between the Fiske steps AVps n = Vn — Ki-i with 
increasing n (cf. Eqs. (fTTj). (O, (fT^ and (j2H)). For large wave numbers, this effect is 
of the same order of magnitude as the reduction of the spacing between Fiske steps due 
to material dispersioriF^. Material dispersion arises from the frequency dependence of the 
complex conductivity and surface impedance of the superconducting elect rodes^^^. It leads 
to a frequency dependence of A/, which becomes significant at frequencies on the order of the 
superconducting energy gap frequency /gap = A//i (A is the superconducting energy gap). 
The gap frequency of the Nb-A10a;-Nb junctions is /gap — 650 GHz (at 4.2 K) and that of 
Al-A10a;-Al is /gap — 85 GHz (at 0.3 K). The frequency dependence of Xl can be written 

AlU) ^ ^ y II - 2/ (. + n.)l (16) 

where /(x) = 1/(1 + e^^^^"^) is the Fermi-Dirac distribution function, ks is Boltzmann's 
constant. 



19 



C. Fiske steps 



Let us consider a long junction of finite length L in the direction of wave propagation. 
Because of refiection at the edges, it behaves like a resonant transmission line and supports 
cavity modes of the electromagnetic field. The theory of such resonances (Fiske resonances) 
in the local case was developed by Kulik^^^. The nonlocal theory for semi-infinite, and 
finite Josephson junctionsiS*^ suggest that for case (c) of Ivanchenko^i, in the presence of 
an external magnetic field Hf., we can use the boundary condition 

= ML) = '^H,. (17) 

Then, the Josephson phase is written as 

(f = ut — kx + ipi{x,t) (18) 

where u = 27tV/^o and k = AnXiHe/^Q. As a first approximation, we consider (fi{x,t) as 
a small perturbation. By inserting this into Eq. (0} we obtain 
d f , 

Aj— / dx'Q{x - x')Lpi^x'ix') - Up^Lfi^tt - ^^p'^(y(fi,t = sm{ujt- kx), (19) 

where we neglect (fi{x, t) in the sine term, and use the assumption concerning the nonlocality 
kernel made at the beginning of Sec. llVyO Furthermore, (pi{x,t) can be expanded in terms 
of the normal modes of the junction 

y^i (x, t) = Im I ^ gj^-' cos k^x I , (20) 

where gn are complex numbers, kn = nir/L. This choice for the x dependence implies the 
boundary conditions v^i,a;(0) = v^i ^.(L) = 0, which is consistent with (fTTj) . The frequency 



ujn is determined from the dispersion relation of Eq. (fTUj) . namely ^uj'^{k) — uj^, uj{k) is 
determined in Sec. IIV Al By inserting the expression (j^Uj) for (fi into Eq. ()19p. and solving 
the set of the equations for gn, we obtain the following expression for the position of the 
Fiske resonances 



K = ^y^^MM^- (21) 

In the local case {Q{x) = 6{x), uj'^{k) = Up + ck'^) we obtain Eq. ©. 

Using Eqs. (HH), (O, and (jH}, in Fig. EUlwe plot the dependence of the first 

Fiske step voltage Vi on the junction width w for AI-AIO2.-AI junctions, and in Fig. IT^ that 
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for Nb-A10j;-Nb junctions. The only fitting parameter for all curves is the Swihart velocity. 
For the Al-AlO^j-Al junctions we found c 8 x 10^ m/s and for Nb-A10a;-Nb junctions 
c 7 X 10^ m/s. The Swihart velocity found using the theory of case^^ (c) and that of 
case^^ (d) differ by about 1%. This small difference indicates that internal nonlocal effects, 
taken into account in case (d), are negligibly small in our junctions. 

In Figs. El and El we compare the measured (symbols) voltage positions Vn of the Fiske 
steps to theoretical curves calculated from Eqs. (fTT|) . (fT^ . (fT^ - (fTB|) and (PT|) . where both 
nonlocality and material dispersion have been taken into account. Solid lines correspond 
to the theory of case^^ (c), and dashed lines to the theory of case^i (d). The data for Al 
junctions is presented in Fig. El and for Nb junctions presented in Fig. The difference 
between the curves according to the two theories decreases with increasing the width of 
the junctions. For the Nb-A10i.-Nb junction of width 4.2 fim the difference cannot be 
distinguished. 

Thus, our experimental data for Fiske steps in both Al-A10a;-Al and Nb-A10a;-Nb junc- 
tions are in good agreement with theory21i24. 



D. Critical magnetic fields and vortex size 

For a junction placed in an external magnetic field the integro-differential equation ^ 
contains an additional Meissner current jm{x, y, z), which is also induced in the absence of the 
junction and depends on the size and the shape of the electrodes (demagnetization effects). If 
the dynamics of the phase inside a small surface region of size ~ ^ L, Aj is disregarded, 
the geometry dependent current jm can formally be taken into account by supplying the 
nonlocal wave equation for the case of infinite length with appropriate (nonlocal) boundary 
conditions. Physically, these boundary conditions reflect edge capacitances and inductances, 
and take into account the stray field in the xy-plane for —w/2 < z < w/2, in addition to the 
stray fields above and below the superconducting leads, which are contained in the nonlocal 
kernel Q{x). These boundary conditions depend on the global geometry of the sample and 
the magnetic screening length in addition to the kernel Q. The magnetic penetration field 
Hci is therefore more difficult to interpret theoretically than the Fiske steps, which follow 
from the kernel Q alone. 

Although in general the solution is difficult to find, in the physically relevant case (c) for 
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Fiske step number 

FIG. 13: Comparison of experimentally obtained Fiske step voltages for Al-AlO^-Al junctions of 
different widths (symbols) with the dispersion curve calculated from nonlocal theories taking into 
account material dispersion (lines). Triangles correspond to the junction of width w = 0.11 /im, 
circles to w = 0.25 nm and rectangles to w = 0.82 /um. The solid line corresponds to the theory of 
case^i (c) and the dashed line to cas©24 (d). 

Al ^ ^1/3, where the Eq. (0} transforms into the local sine-Gordon equation with an effective 
length scale Aj, the expressions for Hd and c = WpXj from the local theory can be used. 
In doing so the correct magnetic screening length, e.g. A^ for w > Xl, and the effective 
vortex size Aj = Aj(l + Xl/wY^"^ from the effective local theory has to be used. Assuming 
that bulk material properties like jc do not depend on the junction width w, it is possible 
to obtain the scaling TelatiomM 

^"^^) - 7^ (22) 

between these quantities, when the width w is changed, which is well satisfied 
experimentally^. Even if the effective local theory provides a good approximation for the 
static case, it might nevertheless fail to describe in the dynamic case some significant prop- 
erties of the junction, such as the possibility of Cherenkov radiation above a critical velocity 
(similarly for the internal nonlocality^) . 
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FIG. 14: Comparison of experimentally obtained Fiske step voltages for Nb-AlO^-Nb junctions of 
different widths (symbols) with the dispersion curve calculated from nonlocal theories taking into 
account material dispersion (lines). Triangles correspond to the junction of width w = 0.3/im, 
circles to w = 1.5 /im and rectangles to w = 4.2 //m. The solid line corresponds to the theory of 
cas^ (c) and the dashed line to case^ (d). 



We will now verify the correct choice of the length scale Aj of a Josephson vortex using a 
variational approach. As the variation parameter we take the vortex size Aj, and minimize 
the Lagrange function of the unperturbed (a = 7 = 0) nonlocal sine-Gordon equation (jlj) 

= j dx - (1 - cosip)^ - J dxdy Lp^{x)ipy{y)Q{x - y) (23) 

(24) 

As trial functions we took the two limiting solutions of Eq. (0]). For w 00, the vortex 
solution of Eq. (j^ is an exponential kink2i 



X 



ip{x, A j) = 4 arctan I exp - 

V Aj 

The other limit is w — > 0, for which the one vortex solution decays algebraically^: 

~ X 
(f{x, Aj) = 2 arctan^3 — h vr. 



(25) 



(26) 
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Here Aj is considered as a free variational parameter corresponding to the characteristic size 
of the vortex. The results for Al-AlO^-Al junctions calculated from the first critical field Hd 
at T = 0.9 K are presented in Fig . El by symbols. The lines correspond to variational results 



for the theories of Refs. 
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and 



2J with the above two ansatzes. Each curve was fitted to the 
experimental data for junction widths in the range 0.7 — 0.8 fim. The prefactors for different 
theories differ by 3 %. Figure ITTl presents the characteristic vortex size Aj as a function of 
junction width for Nb-AlO^-Nb junctions, for T = 4.2 K and Al = 90 nm. In this figure 
the curves are fitted to the experimental data point at width w = 4.2/im. One notes that 
in both cases all four curves overlap with each other. 

The agreement between experimental data (symbols) and theoretical estimations (lines) in 
Figs.^2IHlis rather poor. We suppose this is due to the fact that none of the theories properly 
consider the boundary conditions at the edges of a junction. Deriving appropriate boundary 
conditions for realistic junction geometries in the nonlocal case remains an unsolved problem, 
which requires further work. 

E. Vortex mass 

Finally, we would like to discuss the effect of nonlocal electrodynamics on the dynamical 
mass and quantum behaviour of Josephson vortices. We note that the first experiment with 
a vortex in the quantum regime was presented in Ref. \^ The (nonrelativistic) mass mi;' of 
a Josephson vortex, which appears in the equation of motion 

mpq + ^ = (27) 

for the vortex center of mass coordinate q in an external potential Tiext, is obtained by 
expanding the Lagrange function 

C ^ ^mpv'^ — i^ext + const (28) 

of Eq. in terms of f <^ c. The spacing huo between the energy levels (and, thus, zero 
point fluctuations) is then given by 
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The condition fko^ > ksT is required for dominantly quantum behaviour of the fluxon. In 
the hmit mp 0, the spectrum becomes discrete, provided that the energy scale of the 
potential is fixed. 

The external potential Tiext for a vortex in a Josephson junctions can be, for example, 
induced by an external magnetic field or by presence of a micro-short (micro- resistor) . In 
the first case, the external magnetic field creates a screening current flowing through the 
junction. Provided that the nonlocality of the magnetic screening at the edges can be ne- 
glected or, equivalently, Eq. (fT7|) for the boundary condition is justified, the current spreads 
homogenously over the width of the junction. Hence, the external potential Ticxt is pro- 
portional to the junction width'^^. In the second case, the inhomogeneity leads to a local 
change in the Josephson energy. If the magnitude of the inhomogeneity is constant across 
the junction, the total change in the Josephson energy is proportional to w. As we see, 
in both situations the external potential Tiext scales with w, hence, for characterizing the 
degree of quantum behavior, the use of the specific mass rhp = mp/w for a. vortex is more 
appropriate than using the mass mp- 

For local theory (Aj ^ w ^ A^), the mass of a vortex of form is given by 

8Ejw 

mp = ^T— , (30) 

which vanishes in the limit w —>■ 0, whereas the specific mass rhp = rrip/w is constant. 

Similarly, in the nonlocal case we make use of the ansatz (j){x,t,v) = (f){x,v), where 
X = {x — vt)/Xj{v) and Xj{v) is the characteristic length scale of the vortex, which in 



general depends on the velocity (cf. the Lorentz-factor 1/a/1 — (f/c)^ in the local case), 
and which can be determined from the variational ansatz. As the Lorentz invariance is 
broken in the nonlocal case, we do not know a priori the functional dependence of on v, 
and an additional explicit dependence on v is possible, if the shape of the vortex depends 
on V. From ()23p we obtain 

C 11^07. A 



Ejw 2ul A 
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where 

h = fdx (<^'(x))2; (32) 

dx (1 — cos(^); (33) 

h = - j dxdy (p\xW{y)QM^ - y\)- (34) 

The effective mass in the nonlocal case is obtained by expanding C.{v) into second order. 
Neglecting the explicit dependence on t>, and expanding 




Ajo + T^-, (35) 

one obtains 



Here we again apply the variational approach as we did to obtain the vortex size. Then 
the expression in brackets is the Lagrange equation for the variation parameter Ajo, and is 
therefore equal to zero. Thus the effective mass for the trial function of the form can 
explicitly be written as 

= (37) 

t^pAjo 

Since Aj increases with w, both the effective mass rrif and the specific mass rrif/w (re- 
sponsible for the energy scaling, if ?-^ext ~ w) decrease with decreasing the junction width. 
From Eq. one notes that the inverse of the specific mass scales with the characteristic 
length Ajo (see Figs. Island fTT| . 

Let us emphasize again that the possibility of a vortex behaving as a macroscopic quantum 
particle, when reducing the width w of the junction, not only depends on the scaling of the 
(kinetic) mass rriF, but more exactly the ratio of the kinetic energy to the external potential 
Ticxt; which has to be determined for each specific experimental situation. 



V. CONCLUSION 



The current-voltage characteristics and critical current versus field patterns of narrow 
long Josephson junctions show a strong dependence on the junction width w. This cannot 
be explained using the conventional theory based on the local sine-Gordon model. Our exper- 
imental data are well described by the nonlocal theory originally developed by Ivanchenko^i, 
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and later extended by Lomtev and Kuzovle\«24. According to these models, the electrody- 
namic and static properties of a junction depend on the junction width. Our experimental 
data for Fiske steps observed in both Al-A10a;-Al and Nb-A10a;-Nb long narrow junctions 
are in very good quantitative agreement with these theories. The dependence of the first 
critical field on the width of the junction allows qualitative estimates of the characteristic 
vortex size, which is influenced by nonlocal effects. A proper description of the static effects 
requires a better treatment of the boundary conditions at the edges of a junction. 

We also calculate the specific mass of a vortex in the nonlocal case. According to local 
theory, the specific mass does not depend on the junction width w. By contrast, the nonlocal 
theory predicts that the specific mass of a vortex decreases with decreasing junction width. 
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